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Preliminary results are presented in a step scaling determination of the coefficients in the relativis- 
tic heavy quark action. By matching finite volume, heavy-heavy and heavy-light meson masses, 
we attempt to determine the four parameters (m, £, eg and ce) in the on-shell-improved, heavy 
quark action. In this report we carry out one step in this program by matching two physically 
equivalent systems. The first is a fully relativistic calculation using a 24 3 x48, 1/a = 5.4 GeV 
lattice with both the heavy and light quarks treated as domain wall fermions. The second cal- 
culation uses at 16 3 x 32, 1 /a — 3.6 GeV lattice, a domain- wall light quark and a heavy quark 
computed with the relativistic heavy quark action. The four parameters in this heavy quark ef- 
fective action are then varied to reproduce the mass spectrum from the first calculation. These 
calculations are carried out in the quenched approximation for a heavy quark mass approximately 
that of the charmed quark. 
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1. Introduction 

Heavy quark physics plays an important role in determining the basic parameters of the Stan- 
dard Model [?, ?] and lattice QCD provides a first principles method for determining these pa- 
rameters. However, to treat a charm or bottom quark using presently accessible lattice spacings 
requires the use of an effective theory, one which is not accurate for energy scales on the order of 
the charm or bottom mass but which will accurately describe the physics of charmed or bottom 
states at the energy of scale of interest, the Aqcd scale of non-perturbative QCD. Here we will 
study the relativistic heavy quark effective theory of the Fermilab [?] and Tsukuba [?] groups. 

This relativistic heavy quark action can be written down as: 

S = ^^(x)[m () + 7oA) + C7-5-r,D2-r,£D ; - 2 

x i 

- Y^^OoiFoi + J^^CBOijFij + ^iDo^ijoojWix) (1.1) 

(' hi 

where 1 

Dw{x) = -[V^x)v{x + fL)-Vl{x-p)\r{x-pL)\ (1.2) 

£>JVW = ^[U^x)xi/(x + iX) + ul(x-iX)xi/(x-iX)-2xi/(x)] (1.3) 

Ffi,vY(x) = <T L sJ[U st i(x)Uj v (x + sfl) 

.v,.v'=±l 

xU- s ^(x + sp. +s'v)[/_ji v (i + s'v) — h.c. ]y(*) (1-4) 

As discussed in Ref. [?], this action can be used to compute amplitudes which involve heavy 
quarks carrying spatial momenta of O(p) which will be accurate to order F (ma) (pa) 2 . Potentially 
large errors of order (ma) n , where m is the heavy quark mass, can be removed by a proper choice 
of the seven ma-dependent parameters, m, £, r t , r s , ce, cb and If the coefficient functions F(ma) 
are bounded, then these errors vanish uniformly in the a — > limit. 

As pointed out in Ref. [?], the equations of motion can be used to justify setting r t = 1 and 
^ = 0. With this choice all on-shell Greens functions will still take their continuum form with 
errors of order F(ma)(pa) 2 . As explained in Ref. [?], a final field transformation can be made to 
justify setting the parameter r s = £. Such a transformation will not change particle masses but will 
result in on-shell fermion (e.g. nucleon) propagators which do not show the standard continuum 
form. (Note this field transformation is actually performed on the effective continuum action and 
establishes a relation between the low-energy, on-shell Greens functions of the two theories only.) 
In the work presented here we will make the Fermilab choice, r s = £ and determine the remaining 
four parameters, m, ce, cb and £ using non-perturbative methods that involve only the particle 
spectrum. While not discussed here, the effects of this choice for r s could be determined from the 
spinor structure of the nucleon propagator and then removed. 



2. Step-scaling strategy 

We propose to determine these four coefficients in the RHQ action by matching the finite- 
volume, heavy-heavy and heavy-light spectrum with that determined in an accurate, small-ma cal- 
culation, a strategy similar to that employed for the static approximation [?]. By performing a 
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series of such comparisons, see Fig. [j], we can move from an accurate, ma <C 1 calculation in small 
volume, to a final action at a coarser lattice spacing, practical for large-volume calculations. 



V=24 d x48 



V=16 3 x32 



DWF 

ma=0.2 



a _1 =5.4GeV 

RHQ 

ma=0.3 
a _1 =3.6GeV 




RHQ 

ma=0.45 

a" 1 =2.4GeV 



RHQ 

ma=0.3 
a~ 1 =3.6GeV 

RHQ 

ma=0.45 

a _1 =2.4GeV 



RHQ 

ma=0.65 
a" 1 =1.6GeV 



Figure 1: A schematic of the step scaling technique we use to determine the m, cb,ce and £, appropriate for 
charm quark physics, from non-perturbative, O (a) -improved, light quark calculations at 1 ja = 5.4 GeV. 

Figure |l| shows the steps that we have begun to carry out. Using the quenched approximation 
for simplicity, we match physical quantities calculated on the two finest lattices, which have a fixed 
physical volume of (0.9/m) 3 . For the l/a = 5.4 GeV lattice we use domain wall fermions (DWF), 
involving the single parameter, m ?a 1 .08 GeV, chosen to approximate the charm quark mass. We 
then adjust the coefficients of the RHQ action on the l/a = 3.6 lattice until the spectra of the two 
calculations agree. Second, we expand the volume to (1.35/Vn) 3 , keeping all coefficients fixed, 
and then match with a fourth calculation with 3/2 larger lattice spacing. Here, we demonstrate the 
practical implementation of our approach through stage one: matching the DWF 5.4 GeV (fine 
lattice) spectrum calculated that on a 3.6 GeV (coarse lattice) using the RHQ action. 

In the matching step we compare the pseudo-scalar(PS), vector(V), scalar(S), and pseudo- 
vector(PV) meson masses for heavy-heavy(hh) states and PS and V masses for heavy-light(hl) 
states. We also require the equality of the masses m\ and mi in the hh-PS dispersion relation 
E(p 2 )p h s = m\ + {m\/m2)p 2 ■ For both the heavy quark on the fine lattice and the light quark at 
both lattice spacings we use the DWF action, see Ref. [?] for the method used here and further 
references. In each case we used a fifth dimensional extent, L s = 12 and a domain wall "height", 
M5 = 1.5, ensuring that there are no unphysical propagating states for the quark masses used. 

As a starting point for our l/a = 3.6 GeV calculation, we used the one-loop perturbative co- 
efficients for the five-parameter Tsukuba action [?]. These were translated into the four parameters 
of the Fermilab RHQ action by performing the 0(a) tree-level field transformation: 



y/r = [l + 8fD]xi/ F , \jf T = \if F [i-8YD], 8 



c 



2(m +C)' 



(2.1) 



3. Simulation 



We used the Wilson gauge action since its relation between coupling and lattice spacing has 
been carefully studied [?]. However, as a final check, we examined the static quark potential on 
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our two ensembles. We first determined the lattice spacing directly from the two standards scales, 
ro and rj, obtaining results distorted by our small-volumes. We also obtained the ratio of lattice 
spacings directly from the following relation: 



V l (n l )=V 2 (n 2 /X)/X+C 



(3.1) 



where A(= a%la\) is the lattice spacing ratio between two lattices. The potential Vi{f), from the 
fine lattice, is fit to a standard phenomenological form and then X determined so that this fit, scaled 
as in Eq. [O] matches that on the coarse lattice, giving the expected ratio a\ja% = 1.51(2). Figure. f| 
shows how accurately these scaled static quark potentials agreed. 



I 1 1 1 1 1 1 

5 10 15 

R 

Figure 2: The squares are the static quark potential from the /3 = 6.638 lattice scaled to match the lattice 
spacing at /3 = 6.351. The circles are the j3 = 6.351 potential. 

We employed a "smeared", Coulomb gauge-fixed, heavy quark wavefunction source using 
orthogonal hydrogen ground and first-excited wavef unctions. Table [I] lists the 22 RHQ parameter 
sets that we have run on the coarse lattices. The first few sets of data were chosen close to the 1- 
loop coefficients described above and the latter ones adjusted as we further explored the parameter 
space. 



Table 1: The details of the parameters we have run for the work reported here. 



Label 


m 


CB 


CE 


c 


Label 


mo 


CB 


CE 


c 


1 


0.00 


1.552 


1.458 


1.013 


12 


0.0328 


1.511 


1.438 


1.036 


2 


0.07 


1.547 


1.424 


1.001 


13 


0.0230 


1.511 


1.438 


1.036 


3 


0.0426 


1.550 


1.438 


1.007 


14 


0.0328 


1.511 


1.538 


1.036 


4 


0.0426 


1.550 


1.438 


1.100 


15 


-0.01 


1.700 


1.574 


1.022 


5 


0.0426 


1.550 


1.438 


0.900 


16 


0.00371 


1.709 


1.577 


1.023 


6 


0.0330 


1.609 


1.538 


1.044 


17 


0.0138 


1.715 


1.579 


1.025 


7 


0.0230 


1.609 


1.438 


1.044 


18 


0.02 


1.719 


1.580 


1.026 


8 


0.0230 


1.609 


1.538 


1.044 


19 


0.03 


1.725 


1.582 


1.027 


9 


0.0426 


1.609 


1.438 


1.044 


20 


0.08 


1.707 


1.576 


1.023 


10 


0.0426 


1.550 


1.438 


1.007 


21 


0.09 


1.713 


1.578 


1.025 


11 


0.0 


1.552 


1.438 


1.013 


22 


0.1 


1.719 


1.580 


1.026 
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4. Analysis and result 



We next use the above parameter sets to determine the dependence of the measured quantities 
on the four RHQ action parameters we wish to determine. We attempt to work within a parameter 
range where we can use the linear relation: 



1 coarse 



■ A + J ■ X' RH q 



(4.1) 



Here A is a J-dimensional vector and J an d x 4 matrix. The vector X 1 RH q is formed of the four RHQ 
parameters {(XgnQ)' a }i< a <4 = {m' ,c' B ,c' E ,^'} while the (i-dimensional vector Y l coarse is formed 
from the spectral quantities measured on the coarse lattice. For d = 6 { (Y^ oarse ) a } \ < a <6 = { ^ (Mps + 



3M y )^, {M PS -M v ) hh -\ (M PV -M s ) hh -\ {m y /m 2 ) hh \ \{M PS + 3/V/, 
dex i labels the parameter set and varies between 1 and 22. 



\U,i 



PS 



\hh,i 



}. The in- 



Equation 4.1 can be solved directly for A and J if we use five parameter sets. More parameter 
sets can be used if we minimize a weighted sum of the deviations between the left- and right-hand 
sides of Eq. fO| . We can then use the resulting values for J and A, to solve for the parameters Xrhq 
that yield meson masses as close as possible to those on the fine lattice (denoted as Fy,„ e with error 
Cfine) by minimizing the following quantity with respect to Xrhq- 

(J ■ Xrhq + A — Yfi n e)n \ 



Xfine 



(I 

I 

n=\ 



a 1 

fine,n 



(4.2) 



We divided our data into three categories to explore the sensitivity of the resulting coefficients 
to our choice of masses being matched: A. heavy-heavy system only, including spin-orbit splitting; 
B. heavy-heavy (no SO splitting) and heavy-light systems; and C. all of the above measurements. 
The first five rows of Table ^| demonstrate the stability of our results for the least constrained 
choice of set A. As we incrementally removed the parameter sets that contribute the most to the 
chi-squared in Eq. \2, the matching RHQ parameters remain consistent. As Table § suggests, 



Table 2: A check on the stability of the analysis using the example of data set A and results for B and C. 
Note, the errors on niQ should be viewed as 20% of the w 0.3 renormalized mass. 



data sets 


X 2 /d.o.f. 


m 


cb 


ce 


c 


A(all) 


512/22 


0.035(61) 


1.725(10) 


1.3(5) 


1.036(17) 


(drop 15) 


123/21 


0.034(60) 


1.726(10) 


1.3(5) 


1.036(17) 


(drop 4,15) 


57/20 


0.027(66) 


1.71(9) 


1.3(4) 


1.039(18) 


(drop 3, 4, 15) 


41/19 


0.026(67) 


1.71(9) 


1.3(4) 


1.039(18) 


(drop 2, 3, 4, 15) 


29/18 


0.026(65) 


1.70(9) 


1.3(4) 


1.040(18) 


B 


128/21 


0.03(21) 


1.75(18) 


1.3(17) 


1.04(4) 


C 


42/19 


0.011(61) 


1.73(1) 


1.2(4) 


1.043(17) 



we have very consistent results for different choices of parameter sets and meson masses to be 
matched. The results from parameter set B, which excludes the spin-orbit splitting, gives ce and 
mo with enormous errors. On the other hand, the matching coefficients determined from sets A and 
C look more promising. 
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5. Outlook and Conclusion 

It appeals practical to determine mo, eg, ce and £ from finite-volume, non-perturbative match- 
ing. More statistics and better parameter coverage should reduce the matching errors to a few 
percent. Carrying out the second matching step in Fig. [T] will permit us to perform a RHQ, charm 
physics calculation with \/a = 2.4 GeV and a (2fm) 3 volume. While we have reported quenched 
results here, it should be emphasized that this approach can be extended to full QCD without exces- 
sive computational cost. Since the RHQ parameters being determined are short-distance quantities, 
they depend on both the lattice volume and sea quark masses only through 0(a 2 ) errors. Hence the 
dynamical quark masses need not be physical but must only obey m sea <C I /a. Thus, we need only 
require that Nf = 3 and that m sea /AQCD are equal for each pair of systems being matched. 
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